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Abstract This paper develops a deterministic model of quantum mechanics 
as an accumulation-and-threshold process. The model arises from an analogy 
with signal processing in wireless communications. Complex wavefunctions 
are interpreted as expressing the amplitude and phase information of a mod- 
ulated carrier wave. Particle transmission events are modeled as the outcome 
of a process of signal accumulation that occurs in an extra (non-spacetime) 
dimension. 

Besides giving a natural interpretation of the wavefunction and the Born 
rule, the model accommodates the collapse of the wave packet and other 
quantum paradoxes such as EPR and the Ahanorov-Bohm effect. The model 
also gives a new perspective on the 'relational' nature of quantum mechanics: 
that is, whether the wave function of a physical system is "real" or simply re- 
flects the observer's partial knowledge of the system. We simulate the model 
for a 2-slit experiment, and indicate possible deviations of the model's pre- 
dictions from conventional quantum mechanics. We also indicate how the 
theory may be extended to a field theory. 
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1 Wavefunction analogy in wireless communications 

Several physical systems are characterized by a process of accumulation (of 
energy, charge, etc.), which leads to an activation event once the accumula- 
tion attains a certain threshold. Examples of such systems include lightning 
and nerve impulse transmission. In many cases the accumulation process is 
described in terms of a continuous held, while attaining the threshold trig- 
gers a discrete event. This simultaneous presence of discrete and continuous 
aspects is reminiscent of quantum mechanics. 

Signal acquisition in wireless digital communications also follows this 
same general pattern. Consider a mobile receiver moving randomly within 
a region in which a modulated carrier wave is broadcast. The carrier wave 
is modulated both in amplitude and phase. In order to detect the broad- 
casted signal, the receiver accumulates its received signal until a detection 
threshold is reached. We shall construct a mathematical model of a system, 
in which that the location where detection occurs obeys a probability distri- 
bution reminiscent of the quantum wavefunction. In our model, the wireless 
signal has the following characteristics: 

— The carrier frequency is w, so that the signal has the general mathematical 
form A(r, t) sin(urf + (f>(r,t)). Such a signal is commonly represented by 
its "complex amplitude" A(r,t)e l ^ T ' t \ 

— The transmitted signal (at the transmitter) has constant complex ampli- 
tude over time intervals of length 5, where S » 2-k/uo (S is called the 
"chip width" in digital communications ([1])). The probability distribu- 
tion of complex amplitudes is Gaussian, so that real and imaginary parts 
are independent, identically distributed (i.i.d) standard normal random 
variables with mean and variance 1. 

— The ratio of field amplitude to transmitted signal amplitude (denoted by 
ip(r)) depends on the field location r, but is independent of time. For 
mathematical simplicity, we assume that the ratio assumes one of a finite 
set of complex values {ipi, ■ ■ -ipx}] and that within the (finite) region of 
interest, the sets {r\tp(r) = ipk {k — 1,...K)} all have equal area (see 
Figure 1). 

The receiver has the following characteristics: 

— The receiver consists of an oscillating circuit with natural frequency uj, 
which is driven by the signal field at the receiver's current location. 

— The receiver moves slowly enough so that its field amplitude does not 
change significantly over time intervals of length MS, where M is an 
integer >> 1. 

— The receiver moves in such a way that its position uniformly samples the 
entire region of interest (for instance, by random walk). 

— Our mathematical proof (sec Appendix) requires that the receiver's fields 
over the time intervals (mi, mi + 1)M5 and (m 2 ,m 2 + 1)MS are sta- 
tistically independent whenever mi ^ m 2 . Strictly speaking, a receiver 
moving under random walk will not satisfy this condition: instead, the re- 
ceiver would have to make uniformly-distributed random jumps at times 
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Fig. 1 Wireless communication model 

MS, 2MS, ... A rigorous treatment with random- walk motion would re- 
quire a more careful analysis. 

— The receiver detects the signal when the receiver's amplitude exceeds a 
fixed threshold 0. 

Note this simple model does not include any effects from polarization, propa- 
gation delay or Doppler phase shifting. With the above assumptions, the field 
at receiver position r can be expressed (with the aid of complex amplitudes) 
as follows: 

A(v,t) = M^{v)-v ]t/S] -e iut ] (1) 

where 

ijj(r) = (field amplitude at r) / (signal amplitude at transmitter) 

v\tl$\ may be written as v n = a n + ij3 n , where a„,/3„ are i.i.d. standard 

normal random variables. 
\x] denotes the "ceiling" function, i.e. next largest integer greater than x. 

We now suppose that the trajectory of the receiver is given by the function 
r(t). It follows that the equation for the amplitude x(t) of the driven oscillator 
is: 

x" + uj 2 x = A(r{t),t) (2) 
This equation may be expressed as the real part of the complex equation 

z" +u?z = ^{v{t))-v W5V e i " t (3) 

The solution of (3) which satisfies z(0) = z' (()) = is 
_■ r-t ft 



(4) 

According to our assumptions, the factor e 2luJU in the second integrand oscil- 
lates rapidly compared to the rest of the integrand, which causes the second 
integral to be negligible compared to the first. The model assumptions imply 
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that ip(r(u)) can be treated as constant over time intervals of length MS. 
Using the notation &\ U /(M6)~\ = *P{ r { u )), we have: 

- s r * /<51 

^l^E^/Mr^'^ (5) 

n=l 

The oscillation at time N6 has complex amplitude (—iS/2uj) ■ S(N), where 

N 

S{N) = ^ ln/M] -u n (6) 

n=l 

According to the model assumptions, each \Pj is one of the values {f/'i, • • ■ ~>Pk}- 
Define the random variable n{m) to be the index k corresponding to random 
variable that is 

k(to) ee {k\& m = V'fc} (7) 

Define Nq as the time index at which |5(iV)| first passes a given threshold 
6>: 

N ee nfin{iV| ~ \S{N)\ > 0} (8) 

Our goal is to evaluate the probability distribution of /«(...) corresponding 
to the first passing of the threshold: 

Pv[K(\N e /M]) = k] k = l,...,K, (9) 

This corresponds to the probability distribution of the value of field ip at 
the location of detection. It is a well-known fact in signal processing that 
the rate of accumulation of a random signal is proportional to the signal 
power, which is in turn proportional to the squared signal amplitude ([1]). 
It stands to reason that given a signal that assumes different power levels at 
different times but with equal probabilities, the chance of the accumulated 
signal passing a fixed threshold while at a certain power level should be 
proportional to that power level. This is in fact the case; in the Appendix we 
prove that 

Pr[ K {\N e /M]) =k] oc |^ fc | 2 , k = l,...,K, (10) 
which are exactly the Born probabilities for the spatial wavefunction i[>. 

2 Single quantum detection event model 

In this section, we present a model (based on the model in the previous sec- 
tion) that explains quantum detection probabilities. The model is discretized 
for conceptual clarity and computational tractability; it is fairly straightfor- 
ward to see how the model could be taken to a continuous limit. 

We emphasize that the probability distribution in the previous model 
arose from the outcome of a process. The process involved sampling the 
entire region of potential detection before the actual detection was made. 
This representative sampling was necessary in order for the field strengths to 
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translate into relative probabilities. We want similar characteristics for the 
quantum process. 

In our previous model the process variable is time; this was appropriate 
because we were only concerned about the spatial position of the receiver at 
the moment of detection. However, in quantum mechanics, we are concerned 
about the location of detection events within space-time. It is impossible to 
have a process that unfolds in time that at the same time samples all space- 
time locations before determining the detection location. For this reason, it 
is necessary to introduce a new process variable, so that the process of signal 
accumulation takes place in a non-observable dimension which we will call 
the a-dimension. 

Our wireless communications model had a physical receiver which moves 
within the state space of possible detection locations. Quantum detection 
(say of a particle on a screen) does not appear to have any corresponding 
receiver. We therefore introduce the notion of a detectron, which plays the 
same role as the receiver in our previous model. The accumulation takes place 
as the detectron jumps around and uniformly samples the set of all potential 
detection locations. This "jumping around" takes place in the a-dimension; 
for fixed a, the detectron's space-time location is fixed. We emphasize that 
the detectron is a mathematical construct, and should not be considered as 
a physical particle; we will say more about the physical nature of detectrons 
in Section 5. 

We also postulate a carrier wave that oscillates as a function of a (not as 
a function of time) having the mathematical form sinwa. The frequency lo is 
unknown, and does not correspond to any measurable quantity in space-time. 
The signal has the following characteristics: 

— The signal has constant complex amplitude over a-intervals of length 5, 
where S » 2it/uj. The distribution of complex amplitudes is mean-zero 
Gaussian, with i.i.d. standard normal real and imaginary parts; 

— The signal is multiplied by a complex field amplitude ip(r, t) which is inde- 
pendent of a. For mathematical simplicity, we assume that the amplitude 
assumes one of a finite set of complex values {ipi, ■ • • ipfc}; and that within 
the space-time confines of the detector, the sets {r, t\ip(r, t) — ipk} {k = 
1,...K) all have equal 4- volume. 

The detectron has the following characteristics: 

— Associated with the detectron is an oscillator (which varies sinusoidally 
with a) with natural frequency u, which is driven by the signal field at 
the detectron's current space-time location; 

— The detectron moves in space-time (as a function of a) slowly enough so 
that its field amplitude does not change significantly over a-intervals of 
length MS, where M is an integer >> 1; 

— The detectron moves in such a way that it uniformly samples the space- 
time extent of the detector. 

— The detectron becomes a detection when its oscillator's amplitude exceeds 
a fixed threshold O. 
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We can apply this model to the two-slit setup shown in Figure 2. The de- 
tection moves within the space-time confines of the detection screen. The 
complex field amplitude tp(r,t) corresponds to the conventional Schrddinger 
wavefunction at the screen, which in the ray approximation is given by: 

V>(0, L, z, t) oc dJ 1 e i ^' dl -° J '^ + dq 1 e l{ - k ' d ^'^ (11) 

where k' , u)' are the (observable) wave number and frequency, and d\ — {L\ + 
h 2 y/ 2 + {L 2 2 + {z~h) 2 ) 1/2 , d 2 = {L 2 +h 2 ) 1/2 + {Ll+(z+h) 2 ) 1/2 . We simulated 
this system using MATLAB, with the parameters shown in Table 1 . We 
only considered a single time slice, and restricted to the x — portion of the 
screen. The z locations were discretized into 100 bins; the detectron jumped 
uniformly randomly from bin to bin every M = 400 iteration steps. At each 
iteration, the signal was incremented by v n ■ ip(0, L, z, 0), where v n are i.i.d. 
complex random variables with standard normal real and imaginary parts. 
Each time the detection threshold = 500 was reached, a detection was 
logged and the simulation was restarted. Altogether 100,000 detections were 
logged. Figure 3 shows the detection probability distribution obtained in the 
simulation. Agreement is very close with the theoretical result |?/;(z)| 2 , with 
ip given by (11). Note that in the simulation, distribution peaks are slightly 
lower than theoretical values. If quantum probabilities are indeed the result of 
such an accumulation process, it is possible that measured probability values 
may be lower than the conventional quantum prediction. Unfortunately, it is 
not possible to predict the extent of the lowering from our model, because 
it depends on details of the accumulation process that are not accessible to 
measurement. 



3 Wavefunction formation via accumulation 

The preceding section describes how to obtain quantum- like detection proba- 
bilities given that a certain "broadcast field" is present. However, it provides 
no mechanism for the creation of the broadcast field itself. In this section, 
we show how the broadcast field can be modeled as the result of a pro- 
cess of accumulation that parallels the signal accumulation described above. 
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Fig. 3 Simulation and theory for double-slit experiment with the ray approxima- 
tion 



Table 1 Parameters for double-slit simulation with the ray approximation 



Parameter symbol 


Parameter signifigance (distances in wavelengths) 


Value 


h 


1/2 the distance between slits 


5 


Li 


Distance from source to slit screen 


1 x 10 4 


L 2 


Distance from slit screen to detection screen 


1 x 10 6 


M 


Number from slit screen to detection screen 


400 


-^detect 


Number of detections 


1 x 10 5 


z 


Screen half-width 


1 x 10 6 




Number of bins (discretization) 


100 





Detection threshold 


500 



The well-known path-integral expression for the propagator K(ri, t\\ r 2 , i 2 ) 
is given by ([2]): 

#(!•!, t i; r 2 ,t 2 )= J Dq(t)e lS ^/ h (12) 

Here S[q(t)} is the action, and the notation J Dq(t) denotes an equally- 
weighted summation over all possible paths from (ri,ii) to (r 2 ,i 2 ). This 
integral may be seen as the outcome of an accumulation process. We may 
envision a succession of carrier-wave blips, where each blip corresponds to a 
single path q(t) and makes a differential contribution to the field i/>(r 2 ,i 2 ) 
which is proportional to e lS ^ q ^^ h tp(ri, t\) (as shown in Figure 4). Recall that 
ip(r,t) corresponds to the amplitude and phase of a modulated carrier wave; 
thus e tS ' 9 W/ fi expresses the influence of a source at (ri, t{) on the amplitude 
and phase of the wave at (r 2 , £ 2 ) when a blip passes between them. 

These blips may be associated with detectrons as follows. If (r 2 ,t 2 ) is 
a possible event detection location, then a path q(t) that passes through 
(r 2 ,£ 2 ) can be identified with a detectron location of (r 2 ,i 2 ). We postulate 
that simultaneously with causing an incremental change in the field ip( r 2i ^2), 
the blip also increments the overall complex detection signal amplitude by 
■0(r 2 ,£ 2 ) • v ni as described in the communication model in the previous 
section. In this way, the blips perform a dual mathematical function: they 
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Fig. 4 Single "blip" with field perturbation at (^,£2) 



both build up the field, and furnish the uniform random sampling of possi- 
ble detection sites that is required to obtain Born-rule probabilities (as was 
shown in the previous section) . The process of detection signal buildup and 
detection is shown in Figure 5. 

To see how this works in practice, we focus specifically on the two-slit 
experiment shown in Figure 2 in two space dimensions, using non-relativistic 
electrons of fixed energy fiuj' as particles. We consider the spatial distribution 
of detections at time ti = 0; due to invariance, this distribution will be inde- 
pendent of t 2 - The detection screen corresponds to spatial locations (L 2 ,z). 
We assume the source is configured so that i[)(s+,t) — i()(s-,t) — Ae" % ^ 
where s+ and s_ denote the spatial locations of the two slits ( (0, h) and 
(0, —h), respectively). We also assume that A > is large enough so that the 
signal accumulation process has negligible effect on the size of A. We replace 
t with negative t (since only negative times contribute to detection at t 2 = 0) 
and obtain an expression the Schrodinger kernel for paths that exit through 
the upper slit s + : 

K(s + , -t; r 2 , 0) = k(z - h,t) = (B/t) ■ cxp k)2)) ) (13) 

where B is a constant of proportionality. Similarly, the kernel for paths that 
exit through the lower slit is K(s-, — t; r 2 , 0) = k(z + h,t). The theoretical 
expression for the wavefunction tp(z, 0) is 

pOO poo 

ip(z,0)<x n(z-h,t)e luJ,t dt+ K{z + h,t)e lU3,t dt (14) 
Jo Jo 
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Fig. 5 Detection signal buildup and detection 



which evaluates to 



where 



${z, 0) oc K (-2iy/fUj) + K (-2*V#^) 



Ll + (z±h) 2 ) 
2h/m 



(15) 



(16) 



It was computationally intractable for us to simulate the path integrals that 
give rise to the Schrocdinger kernels themselves. Instead, we assumed the 
kernels and simulated the cumulative effect of different source points. We 
use the following algorithm to accumulate both the fields at different screen 
locations and the overall detection signal: 

Initialize : ip(z) = for all detection locations z; n = 
While | Signal | < 
n = n + 1 
If n divides m 

Change current detectron location z n ; 
Choose random time t (uniformly distributed); 
End If 

^{z n ) = iP{z n ) + n(z n ±h,t)- e*"''; 
v = Normal(0, 1) + i ■ Normal(0, 1); 
Signal=Signal + v ■ ip(z n ); 
End While 

Record detection at location z n 
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Table 2 Parameters for 2-dimensional double slit with non-rclativistic electrons 



Parameter symbol 


Description 


Value (in mks units) 


h 


1/2 the distance between slits 


7.26 x 10~ 10 m 


L 2 


Distance from slits to detection screen 


7.26 x 10~ 9 m 


M 


Number of iterations between detectron jumps 


25 


-^detect 


Number of detections 


3600 


T ■ 
mm 


Minimum value for random time t 


3.37 x 10~ 15 m sec 


T m a x 


Maximum value for random time t 


3.02 x 10" 14 sec 


J 


Electron frequency 


5 x 10 15 Hz 


z 


Screen half-width 


1.45 x 10~ 8 m 




Number of bins (discretization) 


43 





Detection threshold 


1 x 10 5 



Table 2 shows the parameters of the simulation. Figure 6 shows the re- 
sults for 3600 detections. The figure shows counts from 1/2 of the detection 
screen (z < 0), which was discrctized into 43 bins (with the edge of bin 43 
at the center of the screen). The simulation relative counts per bin and the 
accumulated field I^MI 2 are plotted, as well as theoretical detection prob- 
abilities from (15). In general, the theoretical curve lies within error bars of 
the simulation counts; the simulation counts are consistently slightly higher 
than theoretical probabilities near the interference pattern nulls. 



0.08 




Bin Number (Bin 43 is screen center) 

Fig. 6 Accumulated field and detection frequencies for 2-slit experiment 



4 Multiple quantum events 



In the preceding sections we have described how accumulation of the sig- 
nal associated with the moving detectron eventually leads to detection when 



11 



the accumulated signal passes a threshold. This description treats a single 
quantum event in isolation without considering interactions between quan- 
tum events. Accordingly, we now postulate that every space-time event cor- 
responds to a detectron, and all detectron signals at a given a-instant are 
multiplied before accumulation. The situation thus remains as shown in Fig- 
ure 5, except that the jumps in detection signal magnitude are caused by 
products of detectron signals. In fact, the effects of this multiplication are 
already included in the model in Section 2 via the random Gaussian factors 
v n which multiply the held ip(r,t) at the detectron: these fluctuations in the 
signal amplitude correspond to the variation of other detectron signals. 1 

Along with the change in detection signal, at each a-instant the complex 
amplitude of the field at each detectron location is incremented as follows. 
Let (ri, ti), . . . (tm, £m) be the (time-ordered) sequence detectron space-time 
locations for the pre- universe at a. Then the complex wavefunction amplitude 
at (r m ,t m ) is incremented by Y\j <m e tS ^ qjd+1 ^^ n , where gj J+1 (t) is a path 
from (rj,tj) to (rj + i, tj+i). Given that the paths qj 1 ,j 1 +i{t) and qj 2 j 2 +i(t) 
are statistically independent for j\ ^ j 2 , this gives rise to the expression (12) 
for the propagator. 

In summary, each pre-universe shown in Figure 5 gives a single contri- 
bution to the overall detection signal, that accumulates as the pre-universes 
unfold in the a-dimension. Once the detection signal attains a threshold for a 
particular value of a, the space-time universe that we experience is actualized 
as a "snapshot" at that a-instant. 



5 Extension to a field theory 

In the previous sections, we have presented detectrons as moving placehold- 
ers for possible detection sites. This picture presumes a division of the uni- 
verse between fields and detectors. This division of course is not realistic, for 
the detectors themselves are represented by fields in their own right. We may 
eliminate this dichotomy by identifying detectrons as field configurations (in- 
cluding the transmitted particle field and the field representing the atoms of 
the detection screen) that produce contributions to the detection signal. The 
transition to a field theory may then be accomplished by replacing paths 
with field configurations. In other words, the pre-universe at each a-instant 
has a particular field configuration that produces an overall contribution to 
the detection signal. As a increases, the field configurations vary and their 
contributions to the detection signal are accumulated. At the a-instant where 
the threshold is attained, we arrive at the field configuration corresponding 
to the observable universe. 



1 The actual distribution arising from a product of random signals will be log- 
normal rather than normal: we used normal random variables for computational 
simplicity. The results are not sensitive to the particular shape of the distribution. 
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6 Explanations of quantum paradoxes 

Collapse of the wave packet: In our model, the wave function is seen as an 
"actual" field that develops via a process of accumulation. The field is not 
merely a representation of the observer's partial knowledge; but on the other 
hand, the field is not directly observable via physical events in space-time. 
The apparent "collapse" of the wave packet is due to the fact that the physical 
universe is only a single "a-slice" of the entire process. 

EPR and Bell's inequality: Non-local effects pose no problem for this model, 
for the model itself is inherently non-local. An EPR experiment where the 
two spin detectors arc aligned parallel will always detect anti-aligned particles 
because both detections correspond to the same blip and are thus perfectly 
correlated. If the detectors are not parallel, the detection probabilities are 
still determined according to the quantum expression appropriate for that 
configuration. 

Aharanov-Bohm Effect: The Aharanov-Bohm effect shows that fields that 
are localized in a region where a particle can never be detected can still have 
an effect on the motion of the particle. This poses no difficulty to our model, 
because in our model the so-called particle is not an object that travels 
through space-time but rather a correlated series of detection events. 

Identical particle statistics: More work is needed to introduce spin into the 
model. However, in light of our framework it is not surprising that "particles" 
obey special statistics, because "particles" are not separate objects at all. 
What we call a "particle" is simply a series of correlated events. 



7 Comparison With Other Interpretations of Quantum Mechanics 

We briefly compare our interpretation with other alternative interpretations 
of quantum mechanics. 

Everett's "Many-worlds" interpretation^} requires exponential plethoriza- 
tion of space-times. Our model, which embeds space-time within one ad- 
ditional dimension, possesses a much simpler state space. 

Bohm 's quantum mechanics [4] posits that particles such as electrons are able 
to track along with pilot waves. This appears to imply that these particles 
have some sort of inner structure. In our model particles are not "objects" 
at all, so no such complications appear. 

Cramer's transactional quantum mechanics]^} interprets ip* as a wave trav- 
eling backwards in time, but gives no explanation why tpip* should be in- 
terpreted as a probability. Furthermore, transactional quantum mechanics is 
not very clear about the order in which "transactions" are determined. In our 
model, all transactions are determined "simultaneously" (at a = as), and a 
single accumulation process is used to determine all interaction events. 

We also remark that none of these alternative models explains why the wave- 
function is complex, nor why the squared amplitude is interpreted as a prob- 
ability. 
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8 Possible experimental verification 

If true, then our model indicates that the usual formula for a quantum wave- 
function is a statistical approximation, and small deviations from the prob- 
abilities predicted by the wave quation should be expected. In particular, in 
our simulations we consistently found that detection rates near theoretical 
wavefunction nulls were higher that the conventional quantum predictions. 
Unfortunately, the size of these effects would depend on aspects of the process 
that cannot be directly measured. 

9 Conclusions 

Our model presents a radically different picture of reality. Traveling "parti- 
cles" are replaced with series of detection events; as a visual analogy, imagine 
a series of fireflies in a line that flash successively, giving the impression that a 
single firefly is moving along the line. Our model replaces temporal causality 
with atemporal causality; past, present, future are actualized together as the 
result of a process that occurs in a different dimension. Apparent temporal 
causality is due to correlation and not causation. The wavefunction is given 
a physical interpretation as a dynamical field; and the Born rule based on 
the wavefunction is derived as the natural result of a thresholding process in- 
volving this field. The model includes possible differences from conventional 
quantum mechanics. A lowering of peak probabilities in quantum interfer- 
ence patterns compared to the conventional quantum-mechanical prediction 
is a possible result of the model. 

10 Appendix: Mathematical derivation of the Born rule 

In this section, we prove the Born probability rule, 

Pt[ K (\N /M]) = k] cx |^ fc | 2 , k = l,...,K, (17) 

for the wireless communication scenario described in Section 1. We will use 
the notation and definitions of that section. 

In order to investigate the dependence of Pr [k ( \N& /M~\ ) = k] on tp^, for 
each fixed m! > we will investigate the event 

E m ',k = [m' = \N e /M] and k(to') = k] (18) 

conditioned on fixed sequences of (m' — 1) initial ^>'s, corresponding to the 
K m'-i eventg 

F ml ({k' 1 ,...k' m ,_ 1 }) = {k(to) = k' m ,l < m < to', 

where 1 < k' m < K are fixed} (19) 

We shall show that 

Pr\E m ,, k \F' m {{k' x , . . . C_i})] = C(m', {k[„. . . C-i}) ■ IV^I 2 , (20) 
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where C(. . .) is independent of k. The events {F' m {{k\, • • • k' m ,_ 1 })} for fixed 

to' partition the sample space and Pr[F m '({k[, . . . k' m ,_ 1 })] = K 1 ^™ 1 . Fm- 
thermore, the events {£' m / i fe} m , =1 2 partition the event {[k(\No/M~\) = k]}, 
so we obtain 

Pr [k ( \Nq/M~\ ) = k] 

= E E PAE ml!k \F m i{k[,...k' ml _ 1 })]Pv[F ml ({k' 1 ,...k' m ,_ 1 })} (21) 

m' k[,...,k' , 

= I^| 2 E E C(m',k' 1 ,...k' m ,_ 1 )-K 1 - m ' 



(X \lp; 12 



m' k[,...,k' , 



which is the desired result. 

We prove (20) as follows. Conditioned on event ^({/c^, . . . k' m ,_ 1 }) , we 
have 

N 

S(N) = E(VV Fn/M1 ' "n), (JV < to'M) (22) 

71=1 

were the {f„} have i.i.d. standard normal real and imaginary parts (we write 
this as: is n N(0, 1) + iN(0, 1)). It follows that S(N) is a random walk in the 
complex plane with independent (but not identically distributed) steps. We 
also have 

N I 2 

£[|S(A0| 2 ]=2]Tk, (23) 



We shall assume that O » max^ \ipk\, so is much greater than any indi- 
vidual term in S(N). It follows that the distribution of {O^SiN)} for all 
sample paths S can be approximated as a standard Brownian motion B(t), 
where the time variable r is given by 

N I 2 

r(N) ee £ [^-^(AOI 2 ] « 2©- 2 ^ |^ /M1 . (24) 

n=l 

The sample paths comprised in the event |~iV©/M] > m' correspond in the 
Brownian motion picture to sample paths for which \B(t)\ < 1 for all t < 
t((to' - 1)M). For these sample paths, the distribution of _B(t((to' - 1)M)) 
corresponds to the position probability density for a standard Brownian mo- 
tion with absorbing barrier at \z\ = 1. 

Now there is a close connection between Brownian motion and the heat 
equation as follows. Let /3(z, T) be the probability density at time T of a 
Brownian motion with absorbing barrier at \z\ — 1. Then /3(z,T) can be 
found by solving the heat equation with corresponding boundary and initial 
conditions, which in this case are: 

Boundary conditions: f3(z,T) = for \z\ = 1; 
Initial conditions: f)(z, 0) = 5(z), 
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where £(...) is the Dirac delta function. We do not need the complete so- 
lution for f3(z,T) (which can be expressed in terms of the Bessel functions 
{ Jo(a n r/0)}n =1,2,.. .), but we will make use of the following properties: 

a) f3 is radial, so we may write f3(z,T) as f3(r, T) 

b) P(r, T) is C°° for r < and T > 0; 

c) /3 r {l,T) < for all T > 0; 

These properties can be mathematically proven, but are also intuitive conse- 
quences of the physical interpretation of f3(r,T) as an evolving temperature 
distribution within a disk where the boundary is held at zero temperature. In 
light of property c), at time T/3(r, T) can be approximated near the boundary 
\z\ = 1 as 

/3(r,T) = (1 - r)|&(l,T)| +0[(1 - rf]. (25) 
It follows from our identification of {S{N)/0} with B(t(N)) that 

dP[\S((m' - 1) • M)| =r<9 and \N e /M~\ > to'] = A(l-r)+0(l—r) 2 . (26) 



Since \P m = ?/>fc,it follows that the terms &\ n /M] ' v n ~ IV^H-^O, l)+i-N(0, 1)] 
for n = (m' — 1)M + 1 . . . m'M. By rotating in the complex plane we have 
that 



Pr [\S{n)\ < 0,n = (to' - l)M + 1.. .m'M 

\# m , = Vfc and \S((m' - l)M)\ = rO] 



= Pr 



r+\^ k /0\ E 4) 
j=l...J 



< 1, J= 1...M 



(27) 



where ~ AT(0, 1) + i ■ N(0, 1). 

In the case where » max fc \ip k \ and r w 1, we have 



r + \Tpk/0\ E 4 
.7=1.. .J 



Re 



r + |Vfe/0| E 4 
j=i....j 



+ 0([|max|^|/©] 2 ). 



Thus the condition 



(28) 

< 1 can be replaced to a very 



close approximation by the condition Re r + \ipk/&\ £j=i j v 'j 



< 1 (see 
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also Figure 7) and 

Pr[|S(n)| <0,n=(m' - l)M + l...m'M 
\9 m , = V'feand \S((m' - 1)M)| = r0] 



(29) 



Pr 



= Pr 



Re 



Re 



ivw^i E 4 



< 1, J= 1...M 



E 4 
i=i...j 

/©(l-r) 



<G{l-r)/\i) k \,J = l...M 



Note that 4>(x) = 1 for x > J since Re[|i/j|] < 1. Now for sample paths with 





\ 

\ 

\ 

l-r 


fle[z] 


r » 1 


j 

|z| = 1/ 
/ 

/ 
/ 





= 1 



Fig. 7 Depiction of random sum 



\Nq/M~\ > to', the event \Nq/M~\ = m', is the complement of the event 
{IS^n)! < 0, n = (to' — 1)M + 1 . . . m'M}. It follows in summary that 

¥r[E m ,, k \Fm'{{k' 1 ,...k' m ,_ 1 })] 

= Pr [\N e /M] = to' and & m = ^ k \F m ,({k[, . . . 

= Pr [([N e /M] = m'\& m = fa) \F m ,({k[, . . . fc^.J)] • Pr[<f m = ^> fe ] 

-jft-W)) 

• [|S((m' - 1) • M)| = r<9 and ["-2V e /M] > to'] • K~ x . (30) 

Note that (1-^(0(1-^/1^1)) = unless < (l-r) < J ■ \t/j k \/0; and 
since » J ■ \ipk\, the approximation (26) holds on this range. Our integral 
becomes 

«T (wf^^W^l-r)-*. (31) 
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Changing variable to x = (1 — r)/\ipk\i we nave 



«/ (1 - 4>{xG)) ■ |Vfe| 2 • A" _1 Ab • dx (32) 
Jo 

oc \vb k \ 2 . 



Acknowledgements Thanks to Walter Wilcox for many helpful suggestions. 
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